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0. Introduction 

This paper begins a series of papers whose goal is to establish a representation- 
theoretic interpretation of the quantum Knizhnik-Zamolodchikov-Bernard (qKZB) 
equations, and use this interpretation to study solutions of these equations. It was 
motivated by the recent work on the qKZB equations [F,FTVl-2,MV,FV2-5], and 
by the theory of "quantum conformal blocks" that began with the classical paper 
[FR]. 

0.1. The qKZB equations [F] are difference equations with respect to an un- 
known function f(z\, zn, A, t, fi,p) with values in V\ ® ... ® Vjy ® Vjy <8 ■■■ <8> Vj*, 
where Vi are suitable finite dimensional representations of the quantum group U q (Q) 
(q is a simple Lie algebra), Zi,p, r 6 C, and A, /i are weights for g. 

The qKZB equations are a q-deformation of the Knizhnik-Zamolodchikov-Bernard 
(KZB) differential equations, and an elliptic analogue of the quantum Knizhnik- 
Zamolodchikov (qKZ) difference equations, which arc, in turn, generalizations of 
the usual (trigonometric) Knizhnik-Zamolodchikov equations. 

It is proved in [FTV2] (using an integral representation of solutions) that for 
= sl% the monodromy of the qKZB equations is given by the dual qKZB equations, 
which are obtained from the qKZB equations by interchanging (A, r) with (fJ-,p). 
This fact generalizes the monodromy theorems for the KZB and qKZ equations: 
the monodromy of the KZB differential equations is the trigonometric degenera- 
tion of the qKZB equations (which involves dynamical R-matrices without spectral 
parameter) (see e.g. [K]), and the monodromy of the qKZ equations is given by 
elliptic dynamical R-matrices (but there is no difference equation) ([TV1,TV2]; see 
also [FR]). 
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The self-duality of the qKZB equations leads one to expect that they should have 
symmetric solutions Uv 1 ,...,v N (z, X,t, jj,,p), i.e. such that uv 1 ,...,v N {z, X,t, fx,p) = 
u v » v ,(z, fj,,p,X,r), where u* is the dual of u (considered as an endomorphism 
of Vi <g> ... <8> Vn)- Such a solution u (for q = sl 2 ) was constructed in [FV2,FTV2], 
by an explicit integral formula. It is called the universal hypergeometric function. 
This function has many interesting properties, in particular the SL(3, Z)-symmctry 
[FV3-5], where the group SL(3,Z) acts on the lattice 1? generated by the periods 
1 , t, p. A consequence of this symmetry is the qKZB-heat equation [FV2] for the 
function u, which is a q-deformation of the KZB-heat equation [Ber]. 

0.2. A central fact about the KZB and qKZ equations (and one of the main rea- 
sons why they are interesting) is that they are satisfied by conformal blocks. More 
precisely, the KZB equations are satisfied by conformal blocks of the Wess-Zumino- 
Witten conformal field theory on an elliptic curve [Ber], and the qKZ equations 
are satisfied by quantum conformal blocks on the cylinder [FR]. In representation 
theoretic terms, conformal blocks on an elliptic curve are traces of products of inter- 
twining operators for affinc Lie algebras (weighted by an element from the maximal 
torus) [Ber] , and quantum conformal blocks on the cylinder arc highest matrix ele- 
ments of products of intertwining operators for quantum affine algebras [FR]. This 
representation theoretic interpretation of the KZB and qKZ equations is not only 
interesting by itself, but it also allows to prove nontrivial properties of solutions, 
e.g. monodromy theorems (see e.g. [K,FR]). 

The goal of this series is to give a similar interpretation of the qKZB equations. 
In light of the above, the main idea is obvious: one should consider quantum 
conformal blocks on an elliptic curve, or, representation theoretically, traces of 
products of intertwining operators for quantum affine algebras, weighted by an 
element of the maximal torus. It is natural to expect that such traces satisfy a pair 
of dual qKZB equations. This is actually true, and we plan to give a proof of it 
in a subsequent part of the series. However, the details of the proof are relatively 
complicated, and we would like to start with a simpler ("trigonometric") limiting 
case, when t, p — > oo. This limiting case is the main subject of this paper. 

0.3. The structure of this paper is as follows. 

In Section 1 we introduce the main object of the paper - the renormalized uni- 
versal trace function F Vl! ... l v N (X, fi) € [V\ <g> ... <g> V N )[0] <g> {V N <g> ... <g> Vj*)[0], where 
A, p, are weights for q. It is obtained from traces of products of intertwining oper- 
ators for U q (g) weighted by an element of the maximal torus. At the end of the 
section we formulate the main results of the paper - Theorems 1.1-1.5. 

Theorems 1.1 and 1.2 state that the function Fy 1 ,...,v N satisfies two systems of 
difference equations, one with shifts of A, and the other with shifts of fi, which 
go to each other under the transformation A — > fi, /i — > A. In the special case 
= sl n , N = 1, Vi = S mn C n , these systems (as was shown in [EK1]) reduce to 
the trigonometric Macdonald-Ruijscnaars (MR) systems, so we call them the MR 
system and the dual MR system. 

Theorem 1.5 states that the function Fy l! ... l v N (A, fi) is symmetric: Fy 1 .... l v N (X, p) = 
F v * ^y,([j,,X), where * is the permutation of components. It follows from Theo- 
rems 1.1 and 1.2. 

Theorems 1.3 and 1.4 state that the function i*Vi,...,v„ satisfies two additional 
systems of difference equations - the trigonometric degenerations of the qKZB and 
the dual qKZB equations, respectively. 

2 



Theorems 1.1-1.5 are proved in Sections 2-5. 

In Section 6, we study the symmetry of trace functions under q — > q^ 1 , and define 
a modified trace function uy(A, fi), by renormalizcing Fy(A, /i). (This function is 
introduced to connect our paper with the papers [FV2-FV5]; here we define it 
only for N = 1, and plan to define it in general in another paper). Using the 
q — ► g -1 transformation properties and Theorem 1.5, we show that the function uy 
is symmetric. 

In Section 7, we compute the function Fy(A,/z), uy(A, /x) explicitly in the case 
= sk- 
in Section 8, we compute explicitly the trigonometric degeneration of the func- 
tion u from [FV2], in the case N = 1. We show that this function is the same as 
«y(A,//) up to normalization. 

In Section 9, we explain that Macdonald's theory for root systems of type A n -i 
is a special case of the theory developed in this paper, for g = sl n , N = 1, V\ — 

In Section 10, we consider limiting (degenerate) cases of the theory developed in 
this paper. 

0. 4. In subsequent papers of the series, we plan: 

1. To give a representation theoretic-proof of the qKZB heat equation and 
the orthogonality relations for the trigonometric degeneration of the function u 
( [FV2] ) , using the ideas of [EK1,EK2,EK3]. Cherednik's theory of difference Fourier 
transform and Macdonald-Mehta identities for root systems of type A is a special 
case of this theory, corresponding to the situation g = sl n , N = 1, V\ = S mn C n . 

2. To give a representation theoretic derivation of the resonance relations from 
[FV3] in the trigonometric case, using the ideas of [ES] . 

3. To generalize all the results to the case of quantum affine algebras. This in- 
volves a representation theoretic definition of the function u from [FV2] for generic 
values of parameters, for any simple Lie algebra and representations, and a repre- 
sentation theoretic proof of its main properties, such as qKZB and MR equations, 
orthogonality, modular transformations (e.g. the qKZB heat equation). As a spe- 
cial case, this theory should contain Macdonald's theory for affine root systems of 
type A n -\, which was originated in [EK4] but has not been developed from an 
analytic standpoint. In particular, the classical limit (q — > 1) of the modular trans- 
formation of the function u should yield the result of Kirillov [K1,K2] which says 
that the modular transformation S of affine Jack polynomials (which are essentially 
the 1-point functions of the WZW model in genus 1, see [EK4]) is given by a matrix 
of special values Macdonald polynomials at roots of unity. 

4. Specializing this theory to the critical level, to prove that radial parts of the 
central elements of U q {sl n ) at the critical level corresponding to the representation 
S mn C n , are elliptic Ruijsenaars operators (as far as we know, this is known only in 
the trigonometric degeneration). 

Acknowledgments. The first author was partially supported by the NSF grant 
DMS-9700477, and thanks the UNC mathematics department for hospitality. The 
work of the first author was partly done while he was employed by the Clay Math- 
ematics Institute as a CMI Prize Fellow. The second author was supported by the 
NSF grant DMS-9801582, and is grateful to the Harvard Mathematics department 
for hospitality. The authors are grateful to A. Kirillov Jr. for useful suggestions on 
how to improve the paper. 

3 



I. Trace functions for U q (g) 



1.1. The trace functions. 

Let g be a simple Lie algebra over C. Let f) be a Cartan subalgebra of g, and 
oti be simple roots of Q, i = 1, ...,r. Let (ay) be the Cartan matrix of g. Let 
be relatively prime positive integers such that (c^ay) is a symmetric matrix. Let 
ei,fi, hi be the Chevalley generators of g. 

Let t be a complex number which is not purely imaginary, and q = e*. For any 
operator A, we will denote e tA by g -4 . 

Let i7 g (g) be the Drinfeld-Jimbo quantum group corresponding to g. We will 
use the same definition of U q (g) as in [EFK]. Namely U q (g) is a Hopf algebra with 
generators E i} F i7 i~ 1, . . . , r, q h , h E f), with relations: 



q x+y = q x q y , x,y E \) q h E 3 q- h 



1— a» 



EiFj 
1 



q aiih) E t , 

n dihi _ 



q h F 3 q 

-dihi 



-aj(h) p. 



F F 



k=0 
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1—a 



qdi _ q 



EjE? = 0, 



l — dfj 

fe=0 



1 



k 
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where qi = q di and we used notation 



i ^ 3, 



*F 3 F? = 0, i + j. 



[n] q \ 



[k] q \[n-k] q V 



[n] q \ 



[1], • [2], 



in- 



q-q- 



Comultiplication A, antipode S and counit e in U q (g) are given by 
A(Ei) = E t ®q d ^ +1®E U A(Fi) = F t ® 1 + g-* fc < ® F i; A(g h ) = g' 1 ® <?'' 



S{Ei 



-E,q- d ^, S{F i ) = -q d ^F t , S(q h ) = q- 
e(Ei) = e(Fi) = 0, e(q h ) = 1 



Let M M be the Verma module over U q (g) with highest weight /i, w M its highest 
weight vector. Let V be a finite dimensional representation of U q (g), and v G V 
a vector of weight [i v . It is well known that for a generic /x there exists a unique 



intertwining operator $^ : M M — > M M 



V such that = 



w + l.o.t. 



(here l.o.t. denotes lower order terms). It is useful to consider the generating 
function of such operators, e Hom c (M M , ffiM„ <g> F ® V*), defined by = 
J2veB ® w *' wnere the summation is taken over a homogeneous basis £> of V. 

Let Vi, Vjv be finite dimensional representations of U q (o), Vi € vectors of 
weights fj, Vi , such that ^2n Vi = 0. Define the following formal series in 
Vx ® ... ® Vat[0] ® q 2 ( x ^C[[q- 2 ^' ai \ g-2(A,a,)]]. 



(1.1) 



'(A,/x) - ^((^^ ® l^" 1 )...^^) 



It follows from Corollary 3.4, [ES] that this series converges (in a suitable region 
of values of the parameters) to a function of the form q 2 ^ X ' ,1 \f(X,^), where / is 
a rational function in q 2 ^^ and q 2 ^*"*), which is a finite sum of products of 
functions of A and functions of /x. The function (1.1) will be called a trace function. 
Define also the universal trace function, with values in Vi®...®Vjv®V^® ...(giV^*: 

(1.2) *v 1 ...v w (A,m)= E y Vl --' VN (\ri®v* N ®...®vl 

v z eBi 

where Bi are homogeneous bases of Vi. It is easy to see that this function takes 
values in (Vi <S> .. <g> V N )[0] <g> (V£ <g> ... <g> Vi*)[0]. 

Using the generating functions 3>?f , one can express the universal trace function 
as follows: 

(1.3) * yi ... yN (A,/ i ) = Tr((^ Ef=2/i( , i) ® l"- 1 )...*^*), 

where we label the components Vi by i and V™ by *i, and the notation for a 
label k was defined in [F]: when acting on a homogeneous multivcctor, has to 
be replaced with the weight in the k-th component. 

Example 1. (See Section 7) Let g = s?2- In this case let us represent weights 
by complex numbers, so that the unique fundamental weight corresponds to 1. If 
N = 1, and V = Vi is the 3-dimensional representation, then 

<MA,ri = t^i (i + (s 2 - ■»- 2 ) (1 _^ )( CV^-„ ) ) 

(since V[0] is 1-dimensional, we view the function \I/y as a scalar function). This 
example is also computed in [ESI]. 

1.2. The main results. 

It turns out that the trace function \I/yi...v i v( / V/ u ) satisfies some remarkable dif- 
ference equations. These equations are written in terms of so called dynamical 
R-matrices. Below we give a brief introduction to the theory of dynamical R- 
matrices, referring the reader to the expository paper [ESch] for a more detailed 
discussion of them. 

Let V, W be finite dimensional representations of U q (g). Recall from [EV] the 
definitions of the fusion matrix and the exchange matrix. 

The fusion matrix is the operator J wv (/i) : W <8> V — > W <E> V defined by the 
formula 

(1-4) {$™_ llv ®l)$l = <&l wv ^ w ® v ^. 

The exchange matrix Rvw{p) G End(U <£> W) is defined by 

(1.5) Rvwin) = Jywi^ \ 

where 1Z is the universal R-matrix of U q (g). 

We will also use the universal fusion matrix J(A) and the universal exchange 
matrix R(X). They take values in a completion of U q (g) <£> U q (g) and are defined by 
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the condition that they give JywW, RvwW when evaluated in representations 
V,W (cf. also [ABRR],[JKOS]). 

Remark. The fusion matrix describes how to "fuse" together two intertwining 
operators. The exchange matrix describes how to exchange the order of intertwining 
operators. The fusion matrix satisfies the 2-cocycle identity (see below), and is 
sometimes referred to as a "quasi- Hopf twist" . The exchange matrix satisfies the 
quantum dynamical Yang-Baxter equation, and is refereed to as a "dynamical R- 
matrix" . 

Let JJ(A) := J(— A — p), where p is the half-sum of positive roots, and let R(A) = 
R(-X - p). Let Q(A) = m 2 i(l ® S _1 )(J(A)), where m 2 i(a® b) := ba, and S is the 
antipode. It is easy to show that Q(A) is invertible for generic A. 

Define 

(1.6) JJ 1 - Ar (A) = J 1 ' 2 - Ar (A)J 2 ' 3 - Ar (A)..J w - 1 ' JV (A), 

where, for example, JJ 1 ' 2 ---^ stands for (1 ® A„_i)(JJ), where A„_i : U q (g) — ► 
Uqig)™- 1 is the iterated coproduct. We agree that J 1 (A) = 1. Thus JJ 1 -- Ar (A) 
describes how to "fuse" N intertwining operators. 
Let 

(1.7) 5 q (\)=q 2 ^ Jl(l - q- 2 ^ a) ) 

a>0 



be the Weyl denominator. 
Set 



(1-8) Vv 1 ...v w (A,/i) = J 1 - jv (A)- 1 * Vi ...v„(A,aOJ,(A). 

Finally, introduce the renormalized trace function 

(1.9) F Vl ... Vn (X,p) = [Q- 1 ( M )(* Ar )®...®Q- 1 ( M -/ i (* 2 -* Ar ))(* 1 )] Wl ... Vjv (A,-^-p). 

It is convenient to formulate properties of trace functions using this renormalization. 

Example 2. (See Section 7) If q = sl 2 , N = 1, and V = V\ is the 3-dimensional 
representation, then 

_ A g 2(A+ M ) _ q 2\-2 _ q 2^-2 + 1 

F v (\,n) = q 

(This formula is obtained after simplifications from formula (7.20) below when 
m = 1). Note that it is seen from this formula that Fy is symmetric in A and p. 

The following theorems describe the properties of Fv 1 ...v N (A, p). 

For any finite dimensional U q (g)-module W, define the difference operator T>w 
acting on functions of A e f)* with values in (Vi ® ... <g> Vjv)[0] given by the formula 

(1.10) P W = E Tr I^M( R w 1 (A + fc (2 --- W) )...M /y JV (A))T,, 



where T v f(X) = /(A + v), and the component W is labeled by 0. 
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Theorem 1.1. (Macdonald-Ruijsenaars equations) 

(l.H) 'E > wFvi...v N (A, /i) = Xw{q~ 2p ')Fv 1 ...v N (A, (J.), 

where xw{x) — dim W\v\x v is the character ofW, and by T>\j we mean the 

operator T>w acting on F as a function of X, in components V\,...,Vn- 

Theorem 1.1 is proved in Section 2. 

Example 3. If g = sl%, N = 1, and V = V\ is the 3-dimensional representation, 
and W is the 2-dimensional representation, then 

_ (1 _ g 2A-4 )(1 _ g 2A +2) 
W ~ + (l- g 2A-2 )(1 _ g 2A) 

where Tf(X) = /(A + 1). To prove this, it is enough to check that this operator 
is the unique operator of the form a(X)T + b(X)T^ 1 such that (1.11) holds for the 
function Fy given above (note that in our case Xwiq^ 2 ^) = q^ + <Z~ M ). 

Introduce also the dual Macdonald-Ruijsenaars operators D^ Vl acting on func- 
tions of /i e [)* with values in (V^ <g> ... £g> V]*)[0], by the formula 

(1.12) ^ = £ Tr l">] «^ (M + h^-* N -^)...R™ v , {n))T v , 

where V* is considered as a module over U q (g) via the antipode. 
Theorem 1.2. (Dual Macdonald-Ruijsenaars equations) 

(1.13) V^F Vi ... Vn (X,») = X w(q- 2X )F Vl ...v N (\v), 

where T>^ is Vyy acting on F as a function of fi, in components Vfi, V{ . 
Theorem 1.2 is proved in Section 3. 

To formulate the next two results, we need to define some operators acting on 
functions of A and fi with values in (Vi <g> ... ® V N )[0] <g> (V£ ® ... <g> V*)[0]. 
For j = 1, AT define the operators 

(1-14) £>i = (<r 2/i " E X? )*i(«- 2E Xi ® Xi ).j,.i...-j-i. 

where £j is an orthonormal basis of t). Also, define the operators 

Kj = 

(1.15) 

R j+ ij(A + /i^ +2 - JV ')- 1 ...]Rjv J (A)- 1 r j R j i(A + /i^ 2 -^'- 1 ) + /i^' +1 - JV ))...]R jj -_ 1 (A + hV +1 - N >), 

where ^/(A) := /(A + and /iJ- fc acting on a homogeneous multivector has 

to be replaced with the sum of weights of components j, ...,k of this multivector. 
It is easy to check that Dj commute with each other, and its is known that so do 
Kj ([F]). 

Remark. As we have mentioned before, these operators are the trigonometric 
limits of the qKZB operators with spectral parameters, introduced by Felder. 
Analogously, define the operators 

(1-16) D) = (<T 2A -£ *°<)j(q- 2 X Xi ® Xi )j, j+ i...N, 

and 

KJ = R.j-^jifx + /i(* 1 -* J '- 2 ))- 1 ...K,i )H , J -(/i)- 1 r, j x 

(1.17) R. jt * N (n + ftW+ 1 -* JV " 1 ) + h^-* j -^)..M^ j+1 ^ + ht* 1 -*!-^), 

where r»j/(/x) = f(p, + h* 3 ). Like Dj,Kj, the operators Dj ,Kj commute. 
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Theorem 1.3. (the qKZB equations) The function Fy 1 ...v N satisfies the qKZB 
equations: 

(1.18) F Vl ... Vs {\n) = (Kj <g> Dj)F Vl ...v N {\ m)- 

Theorem 1.4. fi/ie c?wa/ qKZB equations) 

(1-19) iV 1 ...v w (A,/i) = (I>y®^y)iV 1 ...v w (A,M)- 

Example 4. Let g = sZ 2 , N = 2, V± = V2 = C 2 with standard basis v + , V-. In 
this case Vi (£> V 2 [0] is 2-dimensional with basis v + (g> w_ , V- <g> w+ , and the action of 
the dynamical R-matrix in this basis is 

k ( a ) = f 

\ q 2 >-~l 

Therefore, if F Vl y 2 (A, /x) is represented by a 2 by 2 matrix with respect to the above 
basis then the qKZB equation corresponding to j = 2 has the form 

(Fu{\,n) Fi 2 (A,/i)\ (V \ 
^21 (A, F 22 (\,h)J \0 q-v) 

(Fu(\ + l,») Fi 2 (A-l,/i)\ 
^2i(A + 1,m) ^(A-l,^)^ 

Here for convenience we took the shift operator from the left side of the equation 
to the right side. 

Remark 1. We should warn the reader that the term "qKZB equations" is nor- 
mally used for equations which contain elliptic dynamical R-matrices with spectral 
parameters, and are difference equations with respect to these parameters z\, zn 
(see e.g.[FTVl]). The equations we consider are a limiting case of the "genuine" 
qKZB equations, when the modular parameter r goes to infinity, and the ratios 
of the spectral parameters Zj/zj+i go to zero, with e - 27rImT << \zj/zj+i\ << 1. 
Namely, the equations considered here are the equations satisfied by the limit (if it 
exists) of a solution of the "genuine" qKZB equations in the described asymptotic 
zone. Throughout this paper, we will abuse terminology and use the term "qKZB 
equations" to refer to this limiting case. 

Remark 2. If N = 1, equations (1.18) and (1.19) are trivial: (1.18) says that 
the t^-component of Fy has zero weight, and (1.19) says that the l / *-component 
of Fy has zero weight. 

Remark 3. It is not hard to show using the quantum dynamical Yang-Baxter 
equation for the dynamical R-matrices that the Macdonald-Ruijsenaars operators 
commute with the qKZB operators. Similarly, the dual Macdonald-Ruijsenaars 
operators commute with the dual qKZB operators. 

Remark 4. Theorems 1.1, 1.3 have the following interpretation. Suppose 
v\,...,Vn are homogeneous vectors in Vi,...,Vjv of weights v\,...,vn, ^ v i = 0- 
Then the function {Fv x ,..,v N (A, fi), v\ ® ... <g> vn) is an common eigenfunction of the 
operators T>^ and Kj with eigenvalues equal to Xw(q~ 2 ^) and 
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Aj(p) = q- 2 ^' v ^+( v i^+ 2 ^i<j^^\ respectively. Thus the trace functions pro- 
vide a solution of the problem of simultaneous diagonalization of the Macdonald- 
Ruijsenaars operators and qKZB operators Kj. 

Remark 5. The problem to deduce equations of type (1.18) for trace functions 
(in the case of affine Lie algebras) was suggested to the first author in 1992 by his 
adviser Igor Frenkel, as a topic of his Ph.D. thesis. However, the first author failed 
to solve this problem at that time, partly because the adequate framework, the 
theory of dynamical R-matrices, was not around yet. 

Theorem 1.5. (the symmetry identities) The function F\r 1 ...v N * s symmetric: 
(1-20) F Vl ,,, VN (X,fi) = Fy* y,([l,,\), 

where F* is the result of interchanging the factors (V\ ® ... ® V/v)[0] and 
(V£®...®1?)[0]. 

Theorem 1.4 is proved in Section 4. Theorem 1.5 is proved in Section 5, using 
Theorems 1.1 and 1.2. Theorem 1.5 and Theorem 1.4 obviously imply Theorem 
1.3. 

Remark 1. Theorem 1.3 can also be derived directly, using the method of 
Frcnkel-Reshetikhin of derivation of the Knizhnik-Zamolodchikov equations. How- 
ever, this derivation is rather long, and we don't give it here. 

Remark 2. In the special case q = sl n , N = 1, V\ = S mn C n , the function F is 
closely related to the kernel of Chcrednik's difference Fourier transform for sl n ([C], 
and the symmetry theorem above (Theorem 1.5) is closely related to Cherednik's 
theorem on the symmetry of the difference Fourier transform [C] . 

Remark 3. In the special case of Remark 2, Theorem 1.5 was proved in [ESI] 
(Theorem 5.6). 

Remark 4. Theorems 1.1-1.5 can be generalized to the case when g is any 
symmetrizable Kac- Moody algebra, and Vj, highest weight modules over U q (g). In 
this case, the functions F make sense as formal power series, and (if g is not 
finite dimensional) the operators T>\y, are infinite difference operators: they 
are infinite sums of terms /(A)T„. However, one can show that these sums make 
sense as operators on power series. We plan to discuss this elsewhere. 

2. The Macdonald-Ruijsenaars (MR) equations 

2.1. Radial parts. 

Let V be a finite dimensional U q (g) module. Then we have the following propo- 
sition. 

Proposition 2.1. (i) For any element X of U q (g) there exists a unique difference 
operator D x acting on V[0]-valued functions, such that 

(2.1) Tr(^Xq 2X ) = D x Tr(^q 2X ). 

(ii) If X is central then Dxy — DyDx for all Y e U q (g). In particular, if X, 
Y are central then DxDy — DyF>x- 

Proof. This is proved in Section 6 of [EK1] (the assumption g = gl n , which is made 
throughout [EK1], is not important for the proof of this result). □ 
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We will call the operator Dx the radial part of X. 

Now recall the Drinfeld-Rcshctikhin construction of central elements of U q (o) 
([D,R]). In this construction, one defines elements Cw corresponding to finite di- 
mensional representations W of U q (g) by the formula 

(2.2) C w = Tr\w(l ® n w )(TZ 21 1Z{l ® q 2p )). 

The map W — > Cw defines a homomorphism of the Grothcndicck ring of the 
category of finite dimensional representations of U q (g) to the center of U q (g). 
Define difference operators Mw '■— P ) c w - 

Proposition 2.2. (i) MwMu = MuMw = Mw®u- 

(ii) M w ^v{\^) = Xw{q 2{M+p) )^v{\n), ^ere X w(x) = dimW[v]x v is 
the character of W . 

Proof. See [EK1] (the proof for gl(n) generalizes tautologically to other Lie alge- 
bras). □ 

This proposition shows that if we can compute A4w explicitly then we will get 
a system of difference equations for $>v 

2.2. The difference equations. 

Let G(A) = q- 2 PQ- 1 (X)S(Q)(\ - h). 

Theorem 2.3. For any V[0]- valued function on \)* , 

(2.3) (M w f)(X)=J2^\wW](G(X + h)R wv (X))f(X + u). 

The proof of Theorem 2.3 occupies Sections 2.3-2.6. 

2.3. The defining property of J. 

Define 

(2.4) J(X) := J{-X- p+ l -{h^ +h^)). 

Lemma 2.4. (The Arnaudon-Buffenoir-Ragoucy-Roche equation, [ABRR], see also 
[3K0S]) Let V, W be finite dimensional U q (g) -modules. Then 

(2.5) ft 2 V A )i JW = JWl E i '® i< (9 2A )i- 

Moreover, this solution is unique among solutions of the form 1 + N , where N has 
only summands with positive weight of the second component. 

Remark. As pointed out in [EV] , this equation is a limiting case of the quantum 
Knizhnik-Zamolodchikov equation of [FR] , obtained when the affine quantum group 
degenerates into the finite dimensional quantum group. 

Proof. We recall ([ABRR], see also [EV], section 9) that J WV (X) satisfies the equa- 
tion 

(2.6) JwvW(q 2{X+p) - E x *h = -R# w q-X*i**itfV+P)--Z *?) 2 J W (A), 
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Using the weight zero property of J, we get 
(2.7) 

Now apply both sides of (2.7) to the subspace W[v] ® V[n]. Using the weight zero 
condition again, we get (on that subspace) 

(2.8) J w (A)(f 2(A+rf )i9f + "(^) 2 - 7^V(9~ 2(A+p) )i(^ + ")i-WA). 
Thus, replacing A with A + |(/x + v), we get 

(2.9) J wv (A + \ (/x + v)) [q-^+P) ) 1 {q v ) 2 = n\} w (q~ 2 ^ ), J wv (A + \ (/x + v)), 

which means that 
(2.10) 

K^q-^+^JwviX+lih^+hW)) 

This implies the first statement. The second statement is straightforward, and also 
follows from [ABRR] . □ 

We will also need the cocycle identity for J. To deduce it, recall the 2-cocycle 
identity for J (see [A,EV]): 

(2.11) J 12 > 3 (A) J 12 (A - h^) = J X ' 23 (A) J 23 (A). 
Thus the cocycle identity for J has the form 

(2.12) J 12 ' 3 (A)J 12 (A + hi^) = J' 1,23 (A) J' 23 (A - h^). 

2.4. The function Zy Let V be a finite-dimensional [7 9 (g)-module. Consider 
the following function with values in V ® V* ® U q (g) (with components labeled as 
1,*1,2): 

(2.13) Z v (\,n)=Tr (^ 01 K 20 q 2X ). 
Lemma 2.5. 

(2.14) Z V (X, fi) = ft 21 (g 2A )iMA, fi). 
Proof. Let us move the R-matrix 1Z 20 around the trace. We get 

Z v (\,n) = Tr a (K 21 K 2 H^ m q ^) = 

(2.15) ft 21 Tr (^°Y A ft 20 ), 

and the Lemma follows after interchangicng q 2X with and moving it around the 
trace. □ 
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Corollary 2.6. 

(2.16) Z v (A,/i) - J 12 (X)^\\ + ^ (2) ,m). 

Proof. Both sides of the equation satisfy (2.14) and have the form {\+\h^ , 
iV, where N has positive degree in component 2. It is easy to show that a solution 
of (2.14) with such property is unique, which implies the lemma. □ 

2.5. The function X v . 

Define the function 

(2.17) AV(A,/x) = Tro^^V^ 03 )- 1 ), 

with values in V ® V* <g> t^(fl) ® t^(fl) (components labeled 1, *1, 2, 3). 
Lemma 2.7. 

(2.18) X v (A,/i) =1l 12 > 3 (q 2X ) 3 X v (\,n)(q- 2X UK 23 )- 1 . 
Proof. We have 

AV(A,/x) - Tr ((^ 03 )- 1 $^ 01 ^ 20 <Z 2A ) = 
^ 13 Tr (^ 01 (^ 03 )- 1 ^ 20 ( Z 2A ) = 
^ 13 ^ 23 Tr (^ 01 ^ 20 (^ 03 )- 1 (? 2A )(^ 23 )- 1 = 

(2.19) ^ 13 ^ 23 (g 2A ) 3 X y (A, M )( g - 2A ) 3 (^ 23 )- 1 . 

The lemma is proved. □ 
Lemma 2.8. 

(2.20) X v {X,n) = J 3 ' 12 (X)J 12 (\ - ^< 3 >)*£>(A + \{h^ - h^)^)J 32 (X)-\ 

Proof. Consider the function X(X) = J 3 ' 12 (A)" 1 X V (\)J 32 {\). (for brevity we 
suppress the dependence on fj, in the notation). By Lemma 2.7 and Lemma 2.4, 
this function satisfies the equation 

(2.21) (g2Aj 3g E(x 4 ®l®x 4 +l®x 4 ®x i )_ X -(A) = X(X)(q 2X ) 3 q^ . 

On the other hand, by the definition of Xy, the function X is of the form 
Zy{\ — i/i' 3 ') + N, where N has negative weights in the third component. As 
before, a solution of (2.21) with such property is unique, so it must coincide with 
its highest term. Thus, N = 0, and by Corollary 2.6 we get the Lemma. □ 

Corollary 2.9. 

Tr (^ 20 (^ 03 )- 1 $^ 01 g 2A ) = 

(2.22) 

(<?- 2A ) 2 ^ 3 ' 12 (A)J 12 (A - h^)^(X + l -{h^ - h^)^)J 32 (\)-\q 2X ) 2 . 

Proof. Straightforward from Lemma 2.8. □ 
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2.6. The difference operators. 

To obtain the action of the operators Mw on traces, we need to apply to both 
sides of (2.22) the operation 77123(1 <B> 1 <£> S), then multiply by q 2p and take the 
trace. 

Let 1Z = Y^ a i®bi, IZ^ 1 = ^2 a i®K- Let u — J2 S(bi)ai be the Drinfeld element. 
Let J{\) = Y.Ck® 4(A), and J -1 (A) = £ <8> d k( X )- Thcn we S et 

Mw*v(A,/i) = 

(2.23) 

Tr 2 ^K a) (A) Cj ® <T 2 ^ 2) (AK(A + ^< 2 >)]*y(A + / l ( 2 ), M )K(A)g 2A 5(4)5(c i ) g 2 ") 2 . 

Let us simplify the expression for £ d' k (X)q 2X S(c' k ), which enters in (2.23). 
Lemma 2.10. 

(2.24) ]T d' k (\)q 2X S(c' k ) = *>P(\)S(u)q 2 \ 

where P(X) := E^A)^ 1 ^). 
Proof. From Lemma 2.4 we get 

J _1 (A)( g - 2A )i =g-^ x '® a! «( g - 2A )iJ- 1 (A)^ 21 . 
Applying rri2i{S <g> 1) on both sides, we get 

(2.25) £4(A) g 2A S(c' fe ) = g Sx?^ d / (A)a . 5(6 . )5(c / )5 2A 5 

which implies the Lemma, since E a iS(bi) = S(u). □ 

Substituting (2.24) into (2.23) and using the cyclic property of the trace, we get 

M w ^v{\v) = 

Y,d^(X)c^ V (X + ^^TT 2 [q^ ^P{\)S{u)q 2X S(cM 2p q- 2X df\\)d J {\+ l - V )] = 
(2.26) 

51 rf| 1 ^(A)c i j*v(A + /j)Tr 2 [(7^ a: *P(A)(7 2A S' _1 (c i )(7 _2A d^(A)d : , (A + i^)5(u)<7 2p ] 
Lemma 2.11. 

^4 1 )(A)® 9 2A 5- 1 (c i ) ( z- 2A 4 2) (A) = 

(2.27) ^(^^W)^" 2 ^ 02 "^ 1 ^ 0255 " 1 ^) 5 " 1 ^) ? 5 ^^)- 
Proof. By Lemma 2.4, we have 

(2.28) ( ( ? 2A )l 1 7^ 23 (A)( (Z - 2A ) 1 = (^23,l)-ljb23 (A)g E(x i )i((x i )2+(x i )s). 

Applying m^S* -1 <g> 1 <g> 1) to both sides, we get the lemma. □ 
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Now we will use the relations 
(2.29) J2 a i (1) S' 1 ^)^ =J2 ak ® ubk > 



(2.30) ^23jl,23 = 

and 

(2.31) ^S^Jdj^A) ® dj 2 } W = S(Q)(A - l -h^) 1 J-\\+ 1 -h^). 

The first relation follows from the coproduct rule for the R-matrix, the second one 
is the functoriality of J, and the third one is obtained if one applies m\2{S® 1 <8> 1) 
to the cocycle relation (2.12). 

Using (2.29)-(2.31), and the fact that (1 <g> S)^' 1 ) = 11, we get from (2.27) 

Y,d?\x)® q 2X s-\ Ci ) q - 2X df\\) = 

^(4 2) (A)a j ) 1 (n- 1 ) 2 g-S(x i ) 2 ((x i ) 1+ (x i ) 2)(5(cfe)d (i) (A)6 . )2 = 

(2.32) q-^ x ^ x ^ + ^\u-%S(m^ - \h)2{J 21 )-\\ + \h {2) )K, 

where h^ 2 > should be replaced by the weight of the second component of the tensor 
product after the action of J^ 1 (i.e. 

J~\\ + \h^)(v ® w) = Ju\^ + \{wt{w) + ® w), 

where wt(w) is the weight of w, and J~ x is the part of J^ 1 which shifts the weight 
of the second component by v). 

Substituting (2.32) into (2.26), and using the identity u~ 1 S(u) — q~ 4p , we obtain 

(2.33) M w *v{\ri=^^\wiv](G(\)R wv (\))* v {\ + v,ii), 
where 

(2.34) G(X) = q~ 2p P(\)S(Q)(\). 

Thus Theorem 2.3 follows from the following Lemma. 

Lemma 2.12. P(X) = Q-*(A + h). Hence G{\) = G{\ + h). 

Proof. The Lemma is obtained by applying mz2i{S~ 1 <g> 1 <g> S^ 1 ) to the cocycle 
identity for J. □ 

Theorem 2.3 is proved. 

2.7. Calculation of G(A). 
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Proposition 2.13. Let G(A) = q-^Q' 1 



(A)5(Q)(A-/i). Then 



(2.35) 



G(A) = 



6 q (X-hY 



The rest of the subsection is the proof of Proposition 2.13. 



Lemma 2.14. 



(2.36) 



A(G(A)) = J(A)(G(A + h™) ® G(A))J- X (A). 



Proof. The element 5(Q)(A) coincides with K'{— A — p), where the element K '(A) 
defined in Section 4.2 of [EV]. By Lemma 27 of [EV], one has (w,K'(X)w*) = 
B\,w(w,w*), where the bilinear form B\,w '■ W ® *W — > C of weight zero is 
defined by the property (1 <g> (,)) o $™_ A , $™* = B\, w (w,w*)Id M> ,- Recall the 
functoriality property of B\,w- 

(2.37) Bu^w ° (Juw <S> </*w*c/) = Bjj o _Bvi/- 

(for suitable A-parameters, which are omitted for brevity). This property implies 
that 



(2.38) A(Q(A)) = {S® 5)(J 21 (A)- 1 )(Q(A) ® Q(A - /i«))J(A - ft* 1 ) - h^)~\ 



which implies the lemma. □ 

Set q = e 7 , Then J(^) is a power series in 7 of the form 1 + 0(7), whose terms 
lie in U q (g) <g> U q (g) <g> F, where F is the field of trigonometric functions. 

Lemma 2.15. Let X(X, 7) = X)m=o^ m (^)T m ^ e a P ower series of weight zero 
satisfying the equation 

(2.39) A(A(A, 7)) = J(A/ 7 )(X(A + 7/^,7) ® *(A, 7))JT 1 (A/ 7 ), 

w/iere Aj e £/ g (g) <8> F 7 swc/i fftai Ao(A) = 1. Then: 

(i) If X = 1 on highest weight vectors of all finite dimensional representations, 
then X = 1. 

(n) In genera/, X = /( {^ 7) > w/iene / = 1 + lh +J 2 f2 + -, ft e F. 

Proof, (i) Xi is a primitive element, so it lies in f) (g> F. Since it acts trivially on all 
highest weight vectors, X\ = 0. Similarly, Xi = for i > 1 by induction, using the 
same argument. 

(ii) Let LOi be fundamental weights, and 771(A) be the eigenvalues of X(A) on the 
highest weight vector of the representation with highest weight Wj. It is easy to see 
that T]i(X + ju)j)r]j(X) = T]j(X + ju>i)rii (A), which implies that there exists g such 
that g(X + -/uJi) = T]i(X)g(X). Let X(A) = X(X)g(X)/g(X+jh). Then X satisfies the 
conditions of (i), so by (i) X = 1. This proves (ii) (we can set /(A) = g(X + jh)). 
□ 
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Corollary 2.16. G(A) = j^z^) f or a suitable function f. 

Proof. Follows from Lemmas 2.14 and 2.15. □ 

To conclude the proof of Proposition 2.13, it remains to show that in Corollary 
2.16, one can take /(A) = S q (X). 

To do this, we apply Theorem 2.3 and part (ii) of Proposition 2.2 in the case 

V = C. In this case \Iv(A,/i) = q g ^ , Rwv = 1- So, using Corollary 2.16, we 
obtain 

This equation is obviously satisfied if / = S. Since the validity of this equation for 
all \l completely determines /(A + h)/f(X), we get the lemma. 

2.8. A modification of Theorem 2.3. 

Let <pv(A,/z) = * v (A,/Lt)5,(A). 

Theorem 2.17. 

(2.41) X>*,¥>v(A,/i) =x W (9 2( " +rt )^(A^). 

Proof. Follows from Theorem 2.3, Proposition 2.2, and Proposition 2.13. □ 

2.9. Proof of Theorem 1.1. 

The theorem follows from Theorem 2.17 and the fusion identity 

(2.42) J 1 - Ar (A)- 1 M ^ 1 - Ar (A)J 1 - Ar (A + &«») = M 01 (A + h^ 2 - N '>)...K 0N (A), 
which is easily checked from the definition. 

3. The dual MR equations. 

3.1. The map Lwv*(n)- 

Let IF be a finite dimensional representation of U q (o). Denote by H^ w the 
isotypic component of type M\ in M M (g) W . For generic fi, we have a natural 
isomorphism r\ : W[X — fi] <g> Ma — ► given by u> <£> w — > 

Consider the intertwining operator 

(3.1) Pv<g>v;wKvw(&% ® %) : M M ® VF -» M M ® FF <g> V ® V*. 

Restricting this operator to the isotypic component of M^ +I ,, from the intertwining 
property we obtain 

(3.2) Pv®v,wT^vw{^ ® Idw)\w[u]®M ll+1 , = L wv » 

where L\yv* {p) W ®V* ^ W ®V* is a, uniquely determined linear map. Our 
task now is to find this linear map explicitly. 

3.2. Computation of Lwv- 
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Proposition 3.1. 

(3.3) L wv .(n) = R wv (n + isy\ 
where ti means transposition in the third component. 

Proof. Let us apply both sides of (3.2) to w <g> y. Let B be a homogeneous basis of 
V. Then we get 

(3.4) J2 Iv^PvwKvw^l ® l)^ +v y] ®v* = ^ W« {»)% +v {™ ®V® «*)■ 

Moving and 7?. to the right, and using the notation Lwv* — X) Pi® Qi> we get 
(3-5) WTt+'Vl ® v * = Pvwt&vw)- 1 E E ^T— t («)*M+^® 

veB i veB 

where wt(v) denotes the weight of v. 
Simplifying (3.5), we get 



(3.6) Yl $X +U){v ® w) y ®"'=EE $rv^ 21 r l J(^)( Pi ™m;v) y g ^ 

This implies that on W[v] (g> V, one has 

(3.7) J( M + v) = U- l J 21 {jx + v){L$ v . . 
That is, 

(3.8) L wv .(ji) = {J- 1 { i i + u)'R 21 J 21 {ii + v))^ = R wv {ii + v) t2 . 
□ 

3.3. Proof of Theorem 1.2. 

Let us multiply both sides of (3.2) on the right by q 2X (acting in ® W), on 
the left by x i\v*i\w ; anc j compute the trace in the Verma modules, summing 
over all v. After multiplication by <5 g (A) and using Proposition 3.1 and the fact 
that ipv has zero weight in V, we obtain 

(3.9) xw{q 2X )^v{\ n) = E Tr| WH (i? w (^ + v) t2 ) Vv {\, fi + v). 
This is equivalent to 

(3.10) X w* {q- 2X )vv{\ n) = E Tr| w . ( _, y) (i? w ( M + z/) tlt2 )^ y (A, /i + v). 

1/ 

Rewriting the last equation in terms of Fy(A, /z), we obtain 
(3.H) 

The expression Wwv(f^) tlt2 can be computed from (2.38): 
(3.12) 

Rwv(v) tlt2 = (Q(^)®Q(^-^ (1) ))Kw*y.(M-^ 1) -^ (2) )(Q" 1 (/-t-^ (2) )®Q" 1 (M))- 

This expression and (3.11) imply Theorem 1.2. 
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4. The dual qKZB equations 

In this section, we will prove Theorem 1.4. 
Using formula (3) in [EV], we get 

(4.1) ® = (i ® t^v^w) E(^-^ fe ® w> 

where 

(4.2) ® w k = Ry W (n)(v ® w). 
Therefore, using simplified notations, we have 

(4.3) *™ +hl v)*v = K-^iM^+Hiw*)*™ = n 2l R*^)-^l +hiw .^ 

Let us now take the j-th intertwining operator in (1.3) and move it to the right, 
permuting it with other operators. We have 

N 

* Vl ... VN (\,ri=n j+lj ...ll Nj (q 2X ) j R* jj+1 (»+ ]T h^)- 1 ...R^)- 1 x 

i=j+2 

Tr((b Vl d)^ -1 (b Vi + 1 3> Vn n 2X <b V i \ — 

N 

n j+lj ...n Nj (q 2X ) j R* jj+1 (n+ ]T /i H) )- 1 -i?* JV (M)- 1 r* J -x 

(4.4) 

1 m+E™ 2 ,^<">' " M+EiLi+i^" 5 H-2;ii j+2 *<•«)• " " 9 M+E^^("> j 

(in the last equality we use that = 0). Using the cyclic property of the 

trace, we can now put the j-th operator into the beginning, and move it to the right 
to its original place, thus completing the cycle. This yields 

"i'u-.ivA.//.. = K l . Kl ...Kx J <<r\ ! K j r...Kj , • 

N N N 

(4.5) *vi...v w (A,/i)- 

Multiplying both sides of (4.5) by S q (X), replacing it by — p — p, and multiplying 
by J 1 "' JV (A) -1 and the product of Q _1 -s, we obtain 

F Vl ... VN (\,ri = 3 1 -- N {\)- 1 n j+lj ...n Nj tf x )^ 

[Q- 1 ^*^ ® ... ® Q _1 (/i - fe (*2-*JV))(*i)] x 

JV N N 

E ^ ( * i) )- 1 -R*iv(M)- 1 r,-/R^(M- E E ^ ( " } 

(4.6) [OM^ ® ... ® Q(/x - /i(* 2 ---* JV ))(* 1 )]Fy 1 ...v N (A, /x). 
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Inverting (4.6), we get 

F Vl ... VN (x,^ = 3 1 -" N (x)- 1 n jj - 1 ...n jl {q- 2X ^ 

[Q" 1 (m) ( * jv) ® ... ® - / i (* 2 -* A, ))(* 1 )] x 

W JV TV 

R Ui(»- E ^ <) )- 1 -^-(m- E /i ( * i) )- 1 r*#^(M)...M* j+1 (M- E ^ 

*=J+1 i=2,ijtj i=j+2 

(4.7) [QM^ ® ... Q(m - / l (* 2 ---* JV ))(* 1 )]Fv 1 ...v N (A, /x). 

Using formula (3.12), it is not difficult to check that on zero weight vectors 

[Q _1 (/i) ( * JV) ® ® - /i(* 2 -* JV ))(* 1 )]x 

AT AT AT 

r ;-i,(m- E E fe ( * i) )- 1 r. J -K; JV (/i)...M; j+1 (/ i - E 

i=j+l i=2,ijtj 

(4.8) [OM^ ® ... ® Q(/x - /j(*2-*JV))(*i)] = x-v. 
So it remains to check that 

(4.9) J 1 "^(A)-%- j _ 1 ...^- 1 (g- 2A ) J -^...^T + i 1 .Ji---^(A) = I>Y 

(on zero weight vectors). 

To check (4.9), observe that it can be written in the form 

(4.10) J 1 - JV (A)- 1 ^,i... j -i( g - 2A )^7 + 1 1 ... Arij J 1 - JV (A) = Dj, 

which implies that it is enough to check (4.9) for TV = 2 and N = 3. 

Let Ej be the left hand side of (4.9). Since E X E 2 E Z = D^D%D% = 1 for N = 3, 
it is enough to prove (4.9) for N = 2 and N = 3, j = 1, 3. But the last two cases 
follow from the case N = 2, so it is enough to check only this case. 

If vY = 2, identity (4.9) easily follows from Lemma 2.4. Theorem 1.4 is proved. 

5. The symmetry identity. 

In this section we will prove Theorem 1.5 using Theorems 1.1 and 1.2. 

Let V be a finite dimensional representation of U g (g), and v <E V a homogeneous 
vector. Consider the operator as a linear operator U q (x\-) — > U q (ri-) <g> V, 
identifying the Verma modules with C/ 9 (n_). 

Lemma 5.1. The operator is a rational function of variables q 2 ^- a ^ of the 
form Q v ^/D{p), where $^ is an operator valued polynomial in Xi := q 2 ^' ai \ and 
D([i) is a polynomial in Xi with nonzero free term. 

Proof. This follows from the arguments in the proof of Proposition 2.2 in [ES]. □ 

Corollary 5.2. The matrix elements of the universal trace function ^ can be ex- 
panded in formal series belonging to the space q 2 ( x ^)C[[q^ 2(X ' ai ^]] <8> C[[g 2 ^' Q ^]]. 

Remark. Note that the tensor product in Corollary 5.2 is algebraic, i.e. un- 
completed. 

Let K = g-2(A,, 1 ) C [[ 3 -2(A,a i )]] ® C[[gr 2 ("' a ')]]. Denote by L the space 
(Vi<g> ...<g> V N )[0]®K. 
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Lemma 5.3. For any vector v G V\ ® ... ® Vjv[0], the function (Fv 1 ...v n ,v) belongs 
to L. 

Proof. It is clear from Proposition 2.2 of [ES] and the definition of J that for 
any V, W the function Ivw(X) is a rational function of yi = q~ 2 ^' ai \ which can 
be represented as a ratio of an operator-valued polynomial of yi whose free term 
is invertible, and a scalar polynomial of with nonzero free term. The same 
statement, for the same reason, applies to Q(A). These facts and Corollary 5.2 
imply the Lemma. □ 

Let L a C L be the space of solutions of equations (1.11) whose highest term (as 
of a series in Xi = q~ 2 ( A >"i)) j s V q- 2 i x ^) ^ where v G (Vi <g> ...Va?)[0] is independent of 
ji. By Theorem 1.1 and Lemma 5.3, for any v G (V\ <g> ... <g> V/v)[0] the inner product 
(F Vl ,...,v N ,v) is in L . 

Lemma 5.4. Lq is a complex vector space of dimension d — dim{V\ <g> ... <£> Vjv)[0]. 
Moreover, if B is a basis of {V\ ® ... <g> Vat)[0] then the collection of functions 
(Fv 1 ,...,v N , v ), v £ B, is a basis of L . 

Proof. The solutions (Fy 1 y N ,v),v G B are linearly independent, since, by the 
arguments in the proof of Lemma 5.3, so are their highest terms. So it remains to 
show that the dimension of the space L is not bigger than d. 

Let I be the maximal ideal in C[[xi, ...,av]]. It is clear that for any m > 1 the 
operator T>w — Xw(q~ 2t *) preserves the subspace I m L, so it is enough to check that 
the dimension of the space of solutions is at most d on L/I m L for all m. 

For s G C*, define an automorphism g s of L by g 8 f{x\ 1 ...,x r ) — f{sx\, ...,sx r ). 
Let Vyy = lim^^o 9sD w9^ X ■ It is enough to show that the system V^f = 
Xw{q~ 2>1 )f has no more than d linearly independent solutions. 

The operator Dyy is easily computed. Namely, since the exchange matrices 
converge to the usual R-matrices as A — > oo (see Theorem 50 of [EV] ) , one gets 

V 

Therefore, the operators are diagonal in the monomial basis, and it is obvious 
that the system V^f — Xw(l~ 2 ^)f nas no more than d linearly independent 
solutions. □ 

Corollary 5.5. There exists a unique element M(p) G End((V^ ® ... <g> V]*)[0]) ® 
C[[q- 2 (M.«i)]] such that 

(5.2) F v ,_ v ,(n,X) = (l®M( M ))iV 1 ,...,v w (A,M). 

Proof. This follows from Lemma 5.4 and Theorem 1.2. □ 

Now we will prove Theorem 1.5. Applying Corollary 5.5 twice, we obtain that 
M(/i)M'(A) = 1 for some function M'(X). This implies that M is in fact indepen- 
dent of fi. Looking at the highest terms of F and F*, one finds that M = 1. 

This proves Theorem 1.5. 
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6. The symmetry of the universal trace 
function under q — » q^ 1 , and the function u v . 

6.1. The symmetry identity with q — > 

In this section we will study the behavior of the functions under the trans- 
formation q — > 5 . 

Consider the algebra isomorphism £ : J7 g (fl) — > l7„-i (g) given by = Eiq dihi , 

£(Fi) — q~ dihi Fi, £,(q h ) — q~ h . Using this isomorphism, we will identify U q (g) and 
f/ 9 -i(g) (as algebras), and regard any module over U q (g) as a module over U q -i(g) 
and vice versa. 

We will start by proving the following theorem. 

Theorem 6.1. The function \I'vi,..,W satisfies the equality 
(6.1) ^v(q-\-X,ti)=u(q)- 1 Q(q,X)^ v (q,X, t l), 

where u(q) denotes the Drinfeld element u ofU q (o). 

Remark. To avoid confusion, here and below we explicitly specify for every 
expression whether it is evaluated at q or at q^ 1 . 

Proof. It is obvious that £ is an antiisomorphism of coalgebras. This implies that for 
a finite dimensional representation V and a vector v € V[0] the operator P$^(q) 
is an intertwining operator — > V <g> M M for [/g-i(g), and hence the operator 
(7£ 10 ) _1 <l?^(g) is an intertwining operator — > <g> U for We have 

(TZ w y 1 (q)^(q)v fi — <g> v + so we have 

(6-2) (7e 10 )- 1 (9)^( g )=^(g- 1 ). 

Let us multiply both sides of (6.2) by q 2X and take the trace. 
Using (2.15) and (2.16), we get 

(6.3) Ti((n w )- l {q)^{ q )q 2X ) = u^Qiq, \)9 v {q, A, /x), 

which implies the theorem. □ 

6.2. The function uy(A,/i). 

Now let V be an irreducible J7 g (g)-module, and define the function 

(6.4) u v (X, fj) := 5 q (\)^ v {q-\ -A, - p). 

Corollary 6.2. The function uy is symmetric: uy{\ A*) = My. A). 

Proof. Let ^ be the highest weight of V. By Theorem 6.1, we have 
(6.5) 

u v (X, M ) - 9~^' I/+2 ' , ^g(A)Q(g, A)* V (A, -/x-p) = «- ( " , " +2p) (Q(«/, A )®Q(<Z, n))F v {\, /*)• 



Therefore, the symmetry of iiy follows from Theorem 1.5. □ 
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6.3. The function uy(X,p). 

The function uy has poles. We would like to get rid of them by multiplying the 
function uy by its denominator. 

Following [ES], introduce the function 

(6.6) Sv(n) = J] JJ(g(°.f«+p)-n(a,a)/2 _ ? -(a, M +p)+r l (a,a)/2^ 

Q >0 n=l 

where for a positive root a we define k% := max{n : F[na] 7^ 0}. 
Define the function 

(6.7) u v (X, p) := S v , (-A - p)5 v (-p - p)u v {\, p). 

Proposition 6.3. The function uy is symmetric: uy(X,p) = u v ,(p, X), and is a 
product of g~ 2 ( A '^) and a trigonometric polynomial of X and p. In particular, it is 
holomorphic in X,p. 

Proof. The symmetry of uy follows from Corollary 6.2. The fact that uy is holo- 
morphic in p is a consequence of Lemma 5.1. The fact that uy is holomorphic in 
A follows from the symmetry. □ 

Using (6.5), one gets the following expression of uy in terms of Fy. 

(6.8) u v (X, p) = g-("."+ 2 ') (Sv. (-X - p)Q(q, X) ® 5 v (-p - p)Q(q, p))F v (X, p). 

In Section 8 we will show that the function uy for U q (sl2) coincides (up to 
a constant factor) with the trigonometric limit of the universal hypergeometric 
function u introduced in [FV2]. 

7. Calculation of the functions u v and F v for sl 2 

7.1. Calculation of the trace function. 

Recall that J7«j(sZ 2 ) is generated by E 7 F,q h , with relations 

(7.1) q h Eq- h = q 2 E, q h Fq-' 1 = q- 2 F, EF - FE = 9 ~ g , , 

q- q- 1 

and the coproduct is defined by 

(7.2) A(E) = E ® q h + 1 ® E, A(F) — F ®1 + q~ h ® F. 

Weights for U q (sl2) can be identified with complex numbers: we say that a 
vector v in an [^(s^-modulc has weight p if q h v = q^v. In this case we have 
(p, p') = pp'/2. We will write q x for q( X ' a \ Thus the meaning of q x in this section 
is different from the previous sections. 

Recall also that for any number a the q- number [a] q is defined by [a] q = q q Z q -i ■ 
Consider the function ^y(X 1 p), where V is the representation of U q (sl2) with 
highest weight 2m, m e Z + . Since the weight space V[0] is 1-dimensional, this 
function can be considered as a scalar function. 
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Theorem 7.1. The function \Iv(A,/i) is given by the formula 

*v(A,/i) = 



(7.3) 

m 

- q,i n^ (i - <? 2 ^-)) n5= (i - <r 2(A - j) ) ' 

where [n] g \ = [l] q ...[n] q . 

Proof. We fix a generator w of V[0]. Let us compute the intertwining operator 
$ M : M„ -» Af„ ® V. 

Let wp, [3 = m,m— 1, —to, be the basis of U defined by the condition Fwp = 
wp-\ if [3 ^ — m (this basis is unique up to a common scalar). The image of the 
highest weight vector under the operator <J> M has the form 



(7-4) $^ = ^ Cj (/x)^ 

1=0 



where Co = 1. 

Let us compute the coefficients Cj(fi). They are computed from the condition 
that A(E) = E (g> q h + 1 (g> E annihilates the r.h.s. of (7.4). Using the formula 

(7.5) ewp = [to + [3 + l] q [m - (3} q wp + \, P^m, 
we can rewrite this condition in the form 

(7.6) q 2i [i] q [n-i + l] q Ci(n) + [m + i] q [m-i + l] g Ci_i(/i) = 0, i > 1. 
Solving this recurrence relation, we obtain 

(7.7, CiW=( _ 1) v«-> B fe|^n[,- i+ it'. 

Now we need to compute $> fl F k v fl . For this, we need to compute A(F k ). Using 
the q-binomial theorem, we obtain 

(7.8) A(F k ) = (F®l + q- h ®F) k = J2 ( k A q~ lh F k - l ®F l , 



where 
(7.9) 



p 

Now, using the intertwining property of (f^, we obtain 
(7.10) (El*) 9"^ fe -'®^)E c ^^ Fj ^ 

Z=0 ^ ^ 9 2 j=0 
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This double sum reduces to a single sum if we use the following version of the 
q-binomial theorem: 

(7.ii) wj) **-'=n(i -A)-*. 

k>l ^ ' P i=0 

Substituting (7.7) into (7.10) and using (7.11), we obtain (7.3). The Theorem is 
proved. □ 

Corollary 7.2. 



(7.12) 



m r _i_ 71 I 

^ q -i(i-i)/2 {q _ q -y rn + l\ , (A+M) jj (<? A + ,_ (J -A-, )( ^_ (Z - Al -, ) _ 



Z=0 Jy j=l+l 

Proof. The statement is obtained from Theorem 7.1 and formulas (6.4) and (6.7). 
□ 

The function uv is manifestly symmetric in A and /x, as predicted by Proposition 
6.3. 



7.2. A formula for 

Now we would like to compute the function Fy . So it remains to compute the 
value of Q(/x) on the zero-weight subspace of V. This value is given by the following 
Lemma. 

Lemma 7.2. The element Q(/i) acts on the zero weight subspace of V by the 
formula 

™ -2^-23+2 _ -2m 
(7.13) Q( M )| y[0] = q ~ 2m II q -2,-2 3 \ • 

J=l 

Proof. Denote by Q r ,z(/x) the eigenvalue of Q(/u) on the subspace of weight I of the 
representation of U q (sl2) with highest weight r. 

It is clear that Qo,o = 1- The values of Q^i and Qi,-i are easily computed from 
the definition. Namely from the ABRR equation we have 

J ip) = 1 + — 5: q \7 q h r-F®E+..., 

q- 2 v — q- 2 (q~ h ® q h ) 



and therefore 

(7.14) <&,_! = 1, Q M 



g 2ai — q 2 
q~ 2 v - 1 



Now consider the subspace of weight I in the tensor product C 2 ® W^, where W 
has highest weight r. Take the determinant of both sides of (2.38) restricted to this 
subspace. Using the strict triangularity of J, we can ignore the J terms and obtain 



(7.15) Qr+l.lMQr-l.lO*) - Q r ,J-l(/i - l)Qr,i+l(/i + 1) 



g- 2 ^ - g- 2 
g- 2 ^ 1 - 1 ' 
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This implies 



(716) Qr+M _ Qr,i-i(/i - 1) q~ 2 >* - q~ 



Q r , I+ l(/i+l) Qr-l,l(/*) <T 2 "-1 ' 

from which we get 

'r+M ir+l TT )/2 ^ V** 1~ r ~ l ~ X 



(7.17) vr+i,tw = TT 

This yields 



r j +1 (p + 1) f = \ q~ 2 ^ +2j - 1 ' q- 2 » - 1 



(r-l+l)/2 ~_ 2j _ 

(7-18) Qh-i,M = II q -2M-^_l ■ 

In particular, 

m „-2n-2j n -2m m „-2p-2j+2 

(7.19, O^M-n 'r^!. -r-n ' r^Vi 

j=0 y j=l H 

as desired. □ 

7.3. Calculation of F V (X, fi). 

Proposition 7.3. The function Fy(X,fi) is given by 



F v (\,ri=q- X »l[ 

(7.20) 



-2n-2j _ 1 



q-2p-2j+2 _ q -2m ■ 



2my^J(i-i)/2 f -iy [m + l] q \ gj^ 

tC {q M^[m-/] g ! n ;. =1 (i-g- 2 ^))n;-=i(i-9- 2 ^)' 

Proof. Using the definition of Fy, Theorem 7.1 and Lemma 7.2, we get (7.20). □ 

Remark. Note that expression (7.20) is not manifestly symmetric. In fact, the 
separate terms in the sum (7.20) arc not symmetric, and it is only the whole sum 
that has the symmetry A — > \i. 

8. Integral representation of the trace function u v for q = sfo 

In [FV2], G.Felder and the second author, studying the qKZB difference equa- 
tions, defined the universal hypergeometric function u m (\,T, (i,p) (depending on 
a parameter q) with a number of interesting properties. In this section we will 
consider the trigonometric limit of this function, and will show that it coincides, up 
to a constant factor, with the function u V ( m ) defined by (6.8), where V(m) is the 
irreducible representation of U q (sl2) with highest weight 2m. 

By definition, the trigonometric limit of u m (A, r, fi,p) is the leading coefficient 
of the asymptotic expansion of u m (\,T, /j,,p) as the modular parameters T,p tend 
to ioo. We will denote this leading coefficient by u m (X, fi). 
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Sending r, p to ioo in the definition of [FV2] , one obtains the following definition 
of the function u m (X, /i). 

Let q = e* with Re(i) > 0, and let < \A\ < 1. Define a function 



, /• " (TjA~ 2 q x+m q-^jTjA-^ 

mi } " q 7|^=...=|T m ,=i (i-r^d-r^) 



- ( I - I - - /';.! - ) 

q-^l-T^f m ,/ T| . 



!<AL (i - ^^- 2 )(i - ) i=1 W=Tr, • 

It is obvious that this function analytically continues to a rational function in A. 
We will denote this analytic continuation also by I m . 

Definition. 

(8.2) «m(A,/i) = I m (X^,q m )- 

Remark. Note that A = q m does not satisfy the condition \A\ < 1, which is 
why we needed to talk about the analytic continuation. 
The main result of this section is 

Theorem 8.1. Let V(m) be the irreducible representation of U q (sl 2 ) with highest 
weight 2m. Then 

(8.3) U V{m) (A, fl) = g (3™-i)™Mi!( g - q~ 1 ) m U m (\, fl), 

where Uyr m -\ is as in Section 6. 

The rest of the section is the proof of Theorem 8.1. 

Lemma 8.2. 

(8.4) 

<\ _ T/T- lN |2 irp 

3 I A m al 3 _ -m(ro-l)/2 



i|T 1 |=... H T m |=i 1< ^ m (l-T l T-V 2 )(l-r 4 T J -V) 27r\/-— TTj [m],!' 



Proof. Let < |p| < 1. Consider the Macdonald denominator of type A m _\ (see 
[Ml]) 

n(x t xj ,p) 
77 -i . . 

where (a,p)<x> := Il^Lo( 1 _ a P J )- 

The Macdonald constant term identity (see [Ml], p. 20-21) says that the constant 
term of the Laurent series (8.5) (with respect to Xi) is given by 

(8 - 6) ci - (A ^ } = ml n fw^' 
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Setting in this identity p = 0, we get 

(i-ty 



(8.7) c.t.(A , t ) = m!- 



Substituting t = g~ 2 , we obtain the Lemma. □ 
Define the expression 

(8.8) I k . m = q -«M» + 2 m) -^M_ Q < fc < ^ 

[k\ q \ 

and the differential form 

k 

n 

3=1 



_ -pr (TjA~ 2 q x+m - q- x - m ){T 3 A- 2 q^ m - q-»~ m ){\ - TjA~ 2 q 2m ~ 2k ~ 2 ) 
k ' m " 11 (1 - T^ 2 )(l - T 3 A- 2 q 2m - 2k ){\ - T 3 A- 2 q- 2 ) X 



^ ,M<k (! - ^ V 2 )(l - TiT^q 2 ) 2nV=lTj 

Lemma 8.3. One ftas 



/ 

J\T 



l\T j \=A 2 q- 2 (^- k )(l+e) 

(8.10) 

kn -2(m-k) (g A+fc - g- A - fc )^ +fc - q-»- k )(l q- 2 ) f 

" (1 - AV 9 <— fc >)(l - ff- a ("-^) J\T^ q -^ fc - 1,m + 



Proof. Let us perform the integration with respect to Tfc, for fixed Ti, ...,Tfc_i. It 
is obvious that the differential form Ffc jm , as a function of Tfe, has two simple poles 
inside the circle of integration -T^ = A 2 q 2k ~ 2m and = 0. Therefore, the integral 
with respect to Tk is equal to the sum of residues at these two poles. The residue 
at the first pole can be found by a direct computation and equals the first term on 
the right hand side of (8.10). The residue at zero equals to Ik, m by Lemma 8.2. 
Lemma 8.3 is proved. □ 

Now let us prove Theorem 8.1. Let us move the contour of integration in the 
definition of I m from \TA = 1 to \TA = A 2 (l + e), via contours \TA = B, A 2 (l + e) < 
B < 1. On the way, we do not run into any poles, therefore, we have 

(8.ii) = <r AM ~ ,n(m ~ 1) / n m . m \ A=qm 

J\ T] \=A 2 (l+e) 

(in the sense of analytic continuation). So, to prove Theorem 8.1, it is enough to 
compute Ji T .i = A 2 (i+e) ^m,m\A=q m ■ We do it by using the recursive relation given 
in Lemma 8.3. Namely, we have 
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Lemma 8.4. One has 

k 

(8.12) / ^fe,mU=(j m = / Ckjlj,m, 

where 

(8.13) cy-(-l) ) 11 (1 _ g2 ( i+m))(1 _ ^(i-^-i)) « 1 ■ 

Proof. The proof is a straightforward induction in fc using Lemma 8.3. □ 
Substituting fc = m in Lemma 8.4, and using the definition of Ik,mi 

we find the 

following expression for u m : 

u m (X^) = q - m <- 3m -^-^ 1 (q- q - 1 )- m x 

(8.14) 

m r i 71 I m 

g -AM ^ g -<( W (*-D/2 (g _ g -iy N +ik (gA+i _ g _ A _ i)(g „ +i _ g _ M _ i) _ 

Comparing (8.14) with (7.12), we get Theorem 8.1. 

9. Trace functions and Macdonald theory 

In this section, following [EK1] and [FV1], we will connect the results of this 
paper with the Macdonald-Ruijsenaars theory. 

We restrict ourselves to the case of g = sl n , N = 1, and let V be the q-analoguc 
of the representation S mn C n . The zero- weight subspace of this representation is 
1-dimensional, so the function <Jv can be regarded as a scalar function. We will 
denote this scalar function by ^ m (q, A, fi) 

Recall the definition of Macdonald operators [M2,EK1]. They are operators on 
the space of functions /(Ai, A„) which are invariant under simultaneous shifting 
of the variables, Ai — > Ai + c, and have the form 

M r = £ (ft ^T_y ) r„ 

/C{l,...,n}:|/|=r \iel,j<£l ) 

where T/Aj = Aj if j ^ / and TiXj = Xj + 1 if j e /. Here q, t are parameters. We 
will assume that t = q rn+1 , where m is a nonnegative integer. 

It is known [M2] that the operators M r commute. From this it can be de- 
duced that for a generic \x = (/Ji, ...,/i„), Mi = 0) there exists a unique power 
series f m o(q, A, /x) G C[[g A2_Al , 9 A " _A " -1 ]] sucn that the series f m (q,X,fi) := 
q2(\,n-mp) f mQ fa \^fj^ satisfies difference equations 

(9.2) M r f m {q,\,n) = ( £ Q 2Ei ^ +p)i )fm(q,X,^- 

/C{l,...,n}:|/|=r 

Remark. The series f m o is convergent to an analytic (in fact, a trigonometric) 
function. 

The following theorem is contained in [EK1]. 
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Theorem 9.1. ([EK1], Theorem 5) One has 



(9-3) f m (q, A, n) = 7m(q, A) 1 ^ m {q \-A,^), 

where 

m 

(9.4) 7m (g,A) :=l[l[(q X '-^ -q 2t q^- X '). 

i=l l<j 

Remark. The exact statement of Theorem 5 of [EK1], in our conventions, 
is that the function f m (q, A, n)-f m (q, A) is equal to T^M^i^ (q~ 1 )q 2X ), which is 
equivalent to Theorem 9.1. 

Let Dw^q 1 ■, ~ A) denote the difference operator, obtained from the operator T>w 
defined in Section 1 by the transformation q — > q^ 1 and the change of coordinates 
A — ► — A. Let A r C™ denote the q-analog of the r-th fundamental representation of 

Corollary 9.2. 

V K r Cn {q-\ -A) = J,(A)7 m («, A) o M r o 7m (< z , A)- 1 ^)- 1 

Proof. This follows from Theorem 9.1 and Theorem 1.1. □ 

In conclusion of this section we would like to make several important remarks. 

Remark 1. Corollary 9.2 is a degenerate (trigonometric) case of Theorem 5.2 
in [FV1], which says that the elliptic Ruijsenaars operators are transfer matrices 
of the elliptic quantum sl n acting in V[0]. Thus, Theorem 1.1 and Theorem 9.1 
immediately imply the trigonometric case of Theorem 5.2 of [FV1] (i.e. the case 
without spectral parameter). 

Remark 2. Conversely, the trigonometric case of Theorem 5.2 of [FV1] together 
with Theorem 1.1 immediately implies Theorem 9.1 (and many other results of 
[EK1]). This is a "direct" proof of Theorem 9.1, in the sense that it involves (unlike 
the original proof of [EK1]) a direct computation of the radial parts of the central 
elements of U q (g). (Another direct proof of Theorem 9.1 is given in [Mi], where 
the radial part of the central element corresponding to the vector representation is 
computed) . 

Remark 3. The line of argument discussed in Remark 2 can be extended to the 
elliptic case. Namely, combining an elliptic analogue of Theorem 1.1 (for affine Lie 
algebras at the critical level), and Theorem 5.2 of [FV1], one can prove an elliptic 
analogue of Theorem 9.1, which says that the radial parts of the central elements 
of U q (sl n ) at the critical level corresponding to evaluation modules A r C n (z), acting 
on functions with values in V[0], are elliptic Ruisjsenaars operators. This has been 
a conjecture for a number of years (see e.g. [Mi], p. 415). We plan to do this in a 
subsequent paper of this series. 

Remark 4. In many arguments of this paper, Verma modules M M can be re- 
placed with finite dimensional irreducible modules L M with sufficiently large highest 
weight, and one can prove analogs of Theorems 1.1-1.5 in this situation (in the same 
way). In particular, one may set ^/ m (q,X,fi) = Tr(<l^g 2A ), where 
: — > Lfj, <£> V ® V* [0] is the intertwiner with highest coefficient 1 (Such an 
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operator exists iff /i — mp > 0, sec [EK1]). Then one can show analogously to Theo- 
rem 9.1 (see [EK1]) that the function f m {q, A, /j.) := "f m {q 1 X)~ 1 ^f m (q~ 1 , —A, p + mp) 
is the Macdonald polynomial P fl (q 7 t,q 2X ) with highest weight /i (/i is a dominant 
integral weight). In this case, Theorem 1.1 says that Macdonald's polynomials are 
eigenfunctions of Macdonald's operators, Theorem 1.2 gives recursive relations for 
Macdonald's polynomials with respect to the weight (for sl(2) - the usual 3-term 
relation for orthogonal polynomials), and Theorem 1.3 is the Macdonald symme- 
try identity (see [M2]). (This representation theoretic derivation of the symmetry 
identity is somewhat different from the one in [EK2] , where a pictorial argument is 
used.) 

10. Limiting cases. 

In this section we will discuss various degenerations of the function i*Vi,...,Viv (q, A, fi), 
and the corresponding degenerate versions of Theorems 1.1-1.5. The main limiting 
cases we will be interested in are the classical limit, and the rational limit. The 
classical limit corresponds to passing from U q (g) to g in the trace construction; in 
this limit the function F depends rationally on p but trigonometrically of A. This 
limit corresponds to the theory of spherical functions on the Lie group G associ- 
ated with g, which is discussed in [EFK1]. In the rational limit, which corresponds 
to the theory of spherical functions on g rather than G, the function F becomes 
rational in both A and p, restoring the symmetry In this limit, the function F is 
the Baker- Akhiezer function for a multivariable bispectral problem (see [Be]). 

10.1. The classical (KZB) limit. 

Let 

(10-1) F^_ VN (X,p) = limF Vu ...y N (q = e\ A M ). 

We will call this limit the classical limit. The existence of this limit follows from 
Proposition 10.1 below. 

Remark. Here and below we write the dependence of functions on q explicitly, 
since in this section q is allowed to vary. 

Example 1. If g = sl(2), N=l, and V — V\ is the 3-dimensional representation, 
we have 

(10.2) ^ (i-i^). 

The classical limit is obtained when in the situation of Section 1 we take the 
ordinary enveloping algebra U(g) instead of the quantized one U q (g). 

More precisely, let 3>Jf be intertwining operators for U(g) defined as in Section 
1, and set 

(10-3) *° Vu ...,Vsfrn) = T^+^m*.) ® l"- 1 )-*^), 

Also, set 

(10.4) S(X)=e^Y[(l-e-^), 

a>0 

and let Q c {n) be the limit of Q(/x) as q — > 1 (i.e. it is defined as in Section 1 from 
representation theory of U(g)). Then we have 
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Proposition 10.1. 

(10.5) 

^...^(A^H^Q- 1 ^)^)®^ 



Proof. The proof is straightforward. 

Let us now look at the degenerations of the properties of -FVi,...,Vn m the classical 
limit. We start with the analogue of Theorem 1.1. 

First of all, we have the following analogue of Proposition 2.1, which is proved 
analogously to Proposition 2.1. 

Proposition 10.2. (i) For any element X ofU(o) there exists a unique differential 
operator Dx acting on V[0]-valued functions, such that 

(10.6) Tr($pfe A ) = L> x Tr($^e A ). 

(ii) If X is central then Dxy = DyDx for all Y e U(q). In particular, if X , Y 
are central then DxDy = DyDx- 

The operator Dx can be computed explicitly for any element X, but in gen- 
eral the answer is complicated. However, if X is the quadratic Casimir C, the 
answer is easy to write down. Namely, define Dx = S(X)DxS(\)^ 1 . Then we have 
([E],[ES1]): ' 



a>0 



where f a ,e a are root generators such that (e a ,f a ) = 1, and A[, is the Laplacian 
on the Cartan subalgebra associated with the standard invariant form. 
Thus, we have the following classical analogue of Theorem 1.1: 

Theorem 10.3. For any X in the center of U(g), let px be the symmetric poly- 
nomial on h* such that X\m„ — px{p- + p)- Then we have 

(10.8) D x F^_ Vn (X^) = P x(-p)F^_ Vn (X^). 

In particular, 



ct>0 



Formula (10.9) was obtained in [E,ES1], but it can also be derived by taking the 
classical limit in Theorem 1.1. 

Now let us consider the classical analogue of Theorem 1.2. Let 2?^ c denote the 
difference operators defined by formula (1.12) for q = 1 (i.e. are the exchange 

matrices for U(g) with \x replaced by — [i — p). Then we have the following result, 
obtained by passing to the limit in Theorem 1.2. 
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Theorem 10.4. 

(10.10) V^F^.yJX^) = Xw{e X )F^ VN {\^), 

Example 2. In the case of Example 1, Theorems 10.3 and 10.4 have the form 

d 2 1 \ u 2 

^F v (X, f ,) = ^-F v (X,f,), 



and 



d\ 2 2 sinh 2 (A/2) J KV 4 



(where T is the shift by 1 in /x), which is easily checked from (10.2). 

Now consider the classical limit of Theorem 1.3. For this purpose introduce the 
classical dynamical r- matrix r(A), which is the classical limit of the exchange matrix 
R(q, A). This matrix is defined by the formula 

(10.11) R(q = e\ A) = i _ 2r(A)t + 0(t 2 ), 

and equals to 

r(A) = --ft + - ^ cotanh-(A,a)e Q A /„, 

where is the Casimir tensor (see [EV]). Taking the quasiclassical limit in Theorem 
1.3, and using that r(— A) = r 21 (A), we obtain the following result. 

Theorem 10.5. For any j = 1, N, one has 

d 



(10.12) [(^ + ^E )*i + E E ^ ® ^°]^,..^ (A, /i), 

w/iere g^yX(A) = J^X(A) i/X is a tensor-valued function whose j-th component 
has weight v . 

The last equation is the trigonometric limit of the KZB equation, which is why 
the classical limit is called "the KZB limit" . 

Let us now consider the classical limit of Theorem 1.4. Let Kj' c be the difference 
operators defined by formula (1.17) for q = 1 (i.e. K(/x) are exchange matrices for 
U(q) with /x replaced by —\i — p). Then we have the following result, obtained by 
passing to the limit in Theorem 1.4. 

Theorem 10.6. One has 

(10-13) Kj' c F Vi Vn (X, fi) = {e x )jF Vi VN (X,fi). 

Finally, Theorem 1.5 does not have an analogue in the classical limit. In this 
limit, the symmetry between A and \i is destroyed, since F c is a product of e^ x -^ 
with a function that is trigonometric in A but rational in /x. 
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10.2. The rational limit. 

The rational limit is a further degeneration of the classical limit. Namely, let 

(10-14) F^_ Vn {\^)=^F- i _ Vn {\ 1 , I x/ 1 ). 

We will call this limit the rational limit. The existence of this limit and the fact 
that det(F r ) ^ can be deduced from Corollary 3.3 in [ESI]. 

Example 3. If g = sl(2), N=l, and V — V\ is the 3-dimensional representation, 
we have 

(10.15) F^(\^)=e-^' 2 (l + i- 

The degeneration of Theorem 1.1 in this limit is the following theorem. Let D r x 
be the rational limit of Dx, i.e. D X (X) is the leading coefficinet of Dx (7 A) as 
7^0. For instance, 



(10.16) D r c = A h -^ 



2 fa & a 



Theorem 10.7. For any X in the center of U(g), let px be the symmetric poly- 
nomial on f)* such that X\m„ = Px(p + P)- Let p r x be the top degree component of 
px- Then we have 

(10.17) D£X,..,^(A,M) =p x {-n)F^...y N {\n). 
In particular, 

(10.18) (A„ - J2 U^)Fv u ...y N {\p) = (p^)F r Vl _ VN (\^). 

a>0 \ ' ' 

The degeneration of Theorem 1.2 looks as follows: 
Theorem 10.8. Equations (10. 17), (10. 18) are satisfied for the function Fyt v » . 
Example 4. In the situation of Example 3, Theorems 10.7, 10.8 have the form 

( d 2 2\„„ x 



which is easily checked from (10.15). 

Using the asymptotics of F v Vn (\,h) at infinity, similarly to arguments of 
Section 5, one can deduce from Theorems 10.7,10.8 the following analog of Theorem 
1.5 (the symmetry theorem): 
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Theorem 10.9. The function Fy v is symmetric: 



(10.19) F v ^ Vn {\,h)=F^I^,\), 

Thus, the symmetry, lost in the first limit, is restored after taking the second 
limit. 

Remark 1. Another proof of Theorem 10.8 is based on representation of the 
above sequence of two limits as a single limiting procedure, which is symmetric in 
A and fi. Namely, one can show that 



(10.20) Fy i y N (X, n) = ]^a^Fv L ,...,v N (q = e s */ 2 , j, ^), 



after which Theorem 10.9 follows from Theorem 1.5. 

Remark 2. Theorems 10.7, 10.8 show that the function F Vi Vjv (A,/i) is a 
solution of the matrix bispectral problem in several variables (on the bispectral 
problem, see e.g. [DG,G]). The Bakcr-Akhiezer function of the rational Calogcro 
system of type A, which is a known solution of the multidimensional bispectral 
problem ([VSC], see also [Be]), is a special case of Fy VN (\,fi) (N = 1, V = 
Vi = S mn C n ). 

Finally, let us consider the rational limit of Theorems 1.3, 1.4. To formulate 
the analog of Theorem 1.3, introduce the rational limit of the classical dynamical 
r-matrix, r°(A) = lim 7 ^ l r (l^)- It has the form 

(10.21, r » (A) = E ^ 

Taking the rational limit in Theorem 10.4, we get 
Theorem 10.10. For any j = 1, N, one has 



(E^w-E^w) 



(10.22) ^F Vi _ Vn {\^). 

The analogue of Theorem 1.4 is 
Theorem 10.11. Equation (10.22) is satisfied for the function Fyt v » . 

10.3. The qKZ and KZ limits. 

Assume that |g| < 1. The qKZ limit is defined by 

(10.23) F^l (A,/i)= lim q 2 ^F Vl _ VN {\^). 

{\,a.i)—> — oo 

It is easy to check using Theorem 50 of [EV] that 

Fy«l,y N (X,n) = 

[Q- V)^ ® ® Q~V - ^ 2 -* JV) ) ( * 1) ]<(^Ef- 2 ^» ^ N - 1 )-K N ) 

(10.24) [Q~ 1 (v) i * N) ® ••• ® Q~V - /i ( * 2 -* JV) ) ( * 1) ]J 1 - JV (m)*, 
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where (, ) denotes the highest matrix element. (The last expression is an endomor- 
phism of (Vx <8> ... ® Vf")[0], which is regarded as an element of 
(Vi <8> ... ® Vjv)[0] ® (V$ ® ... ®V*)[0]). In particular, this function is indepen- 
dent on A. 

Let us now consider the behavior of the equations given by Theorems 1.1-1.5 in 
the qKZ limit. 

The MR equations given by Theorem 1.1 become trivial. Namely, when (A, a,) — > 
— oo , one has J(A) — * 1, and hence K(A) — > 7? 21 . The matrix 1Z 21 is triangular, so 
only its diagonal part contributes to the trace. Inspection of this diagonal part 
shows that lim/x, a A—y—oo "Dw = Y) t . dimW\u] T v , and the limiting equation is 

(5^dimW[i/] q -^)T v )F qKZ = xw{q~ 2 ^)F qKZ 
v 

(where T„ is the shift of A), which is a trivial consequence of the fact that F qKZ is 
independent on A. 

The dual MR equations given by Theorem 1.2 have a slightly more interesting 
limit. It is easy to see that the only term on each side of (1.11) which survives in 
the limit is the term corresponding to the lowest weight vw of W . Therefore, the 
limiting equation has the form 
(10.25) 

The qKZB equations become the trigonometric limit of the qKZ equations. 
Namely, for all j — 1 , . . . , N we have 

(10.26) [KT] +v ..nji(q- 2 n j K lj ...n j - ltj ® ^]<^ y >) = F vf, Z ,V»- 

If N = 2, these equations are closely related to the Arnaudon-Buffenoir-Ragoucy- 
Roche (ABRR) equation (see Lemma 2.4). In general, they are essentially the 
N-component version of the ABRR equation. 

Remark. The dual qKZB equations do not seem to have a reasonable qKZ 
limit. Also, the symmetry relation (Theorem 1.5) does not hold in the qKZ limit 
since the function depends on ft and not on A. 

The KZ limit is obtained from the qKZ limit as q — > 1, in which case the qKZ 
equations degenerate into the trigonometric KZ equations (the quasiclassical limit 
of the ABRR equation). We leave it to the reader to derive the limiting equations 
in this case. 

We plan to consider these limits in more detail in another paper in the more 
interesting case of afhne Lie algebras and quantum affine algebras. 
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